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I. INTRODUCTION 

Biphoton states generated in Spontaneous Parametric Down-Conversion (SPDC) can be 
highly entangled in continuous variables characterizing signal and idler photons, such as 



angular variables or frequencies 



y, |2j. For SPDC processes with the type-I phase matching 
the degree of entanglement is characterized by values of the Schmidt number K up to 
300 in the case of short pump pulses. Such high entanglement reflects a high degree of 
correlations accumulated in the biphoton state and, in principle such highly-correlated states 
can be useful for goals of quantum information and quantum cryptography. But, of course, 
there is a big gap between finding that the degree of entanglement is very high and its 
application in practice. We think that an important step in the direction of closing this gap 
is finding the Schmidt modes and eigenvalues of the reduced density matrix rather than only 
the Schmidt number K. This task is not quite trivial because Schmidt modes are known 
only for double-Gaussian wave functions 4|, whereas in some cases even a possibility of 
modeling the biphoton wave function by a double- Gaussian function is problematic. We 
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will show below that in such reasonable way of derivation consists in calculating 

first the reduced density matrix (even if approximately) and, then, modeling the latter by a 
double-Gaussian function. Eigenfunctions of such a model double- Gaussian density matrix 
are just the Schmidt modes. By continuing this procedure, and using general connections 
between the density matrix and wave function, we show that the double-Gaussian model of 
the spectral wave function exists, and it's found explicitly. 

The paper is structured in the following way. In the next Section we give a brief overview 
of the relation between the double- Gaussian wave functions and density matrices and deriva- 
tion of their Schmidt modes. In Section 3 we summarize the main features of the spectral 
biphoton wave function for the degenerate coUinear SPDC process with a pulsed pump. In 
Section 4 we describe how the Gaussian modeling can be applied for finding Schmidt modes 
in the case of long pump pulses. In Section 5 we describe Gaussian modeling of the re- 
duced density matrix and derivation of the Schmidt modes in the case of short pump pulses. 
In Section 6 we formulate some generahzations of the derived specific results and suggest 
a single double- Gaussian model of the biphoton spectral wave function valid for arbitrary 
pump-pulse durations. 



II. REDUCED DENSITY MATRIX AND SCHMIDT MODES FOR 
DOUBLE-GAUSSIAN BIPARTITE WAVE FUNCTIONS 



Let Xi and X2 be some continuous variables of particles "1" and "2" and \l/(xi, X2) be 
a bipartite wave function. Lets us consider here only the case of symmetric bipartite wave 
functions for which "^{xi, X2) = '^{x2, Xi). The most general form of the double-Gaussian 
wave function satisfying this condition is given by 

(xi + ^2)^ 



^'(Xi, X2) 



irab 



exp 



2a2 



exp 



(a^i - X2) 
262 



(1) 



where a and b are parameters related to widths of distributions in xi^2 determined by each 
Gaussian factor in Equation ([I]), a > and 6 > 0. Alternatively, the same wave function can 
be presented in a different form, convenient for determining the double- Gaussian Schmidt 
modes {4] 



a 



\E'(xi, X2) = exp < 



2\ (^1 + ^2) + 



2/i 



2(1-/^2 



X\X2 



(2) 
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The parameters a and b are related to /i and a by equations 
a — b 



a + b 



a = A / — and a = — i / , o = — i / 

ao a \ 1 — fi a \ 1 + fi 



(3) 



-1 < /i < 1. 

The density matrix and reduced density matrix, corresponding to the wave function 
^(a^i, X2), are given by 



p(xi, X2; x[, X2) = "^{xi, X2)'^{x[, X2) and 



Pr{Xl^ X-^j 



dX2p{Xi, X2] x\, X2). 



(4) 



After integration, the reduced density matrix can be presented in the same forms as described 
above for the wave function, ([T]) and ^ 



Pr{Xl, Xi^ 



1 2 



{xi+x'^f 




(xi-x;)2" 


exp 


252 


262 . 
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+ P /~2 




1 


— — exp < a 




-/i t 



1 + /^ / 2 , ,2. 



[X^ + Xi ) + 



2J1 



(5) 



The pairs of parameters a, 6 and 5, /I are related to each other by the same formulas as a, b 
and a, /i (j3]): 

^2 /l + M r V2 



P 




and a 



- P 



(6) 



a + 6 V a6 « 

Moreover, with a simple algebra we find a series of formulas relating the parameters char- 
acterizing the density matrix with the parameters characterizing the wave function 



a = a, p = p , a 



62 



ab 



+ b^' 



(7) 



and 



2 ' 



a 



a — b 



(8) 



The Schmidt modes {■?/'„ (xi ), Xn (2^2)} are defined as eigenfunctions of the reduced density 
matrices Pr{xi, x[) and Pr{x2, X2). For symmetric wave functions X2) two functions 

in the Schmidt-mode pairs are identical, Xn{x2) = i'n{x2). In a general case, to find Schmidt 
modes one has to solve the integral equation 



dx[pr{Xi, x[)iJn{x[) = Xn1pn{xi), 



(9) 
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where A„ are eigenvalues of the reduced density matrix. The Schmidt modes diagonahze the 
reduced density matrix and, hence, the latter can be presented as the sum of products of 
Schmidt modes 

oo 

Pr{Xi, X[) = ^ \ni'n{Xl)i^nix[) . (10) 
n=0 

According to the Schmidt theorem, the wave function "^{xi, X2) also can be expanded in a 
series of products of the Schmidt modes 



(Xl, X2) = ^ a/A^ ^niXi)^IJnix2). 



(11) 



n=0 



For the double- Gaussian wave functions analytical expressions for the Schmidt modes are 
known and they follow from the formula ^ 



exp < a 



-{xi + X2] 



2/i 



2(l-/i2)^ (l-/i^ 
^/n ^/T^^^(^''^Pn{aXl)^Pniax2), 



X1X2 



(12) 



n=0 



where are the oscillator wave functions ipn{x) = {2"'n\y/7T)^^^'^e~^^ Hn{x) , Hn{x) are the 
Hermite polynomials. The Schmidt modes and eigenvalues of the reduced density matrix 
are easily found now from Equations ([I]), ([2]), ([3]), and ( |T2l) 



and 



i^n{xi^2) = y/a ipn{axi^2) 



ab 



\ab 



2v ab f a — b 



a + b \a + b 



-1)"1 



2b a-b 



n/2 



a + b \a + b 



2n 



A„ = 4a6 



2b a-b 



{a - bf 



(13) 



(14) 



The Schmidt number, characterizing the degree of entanglement in the state ^^(xi, X2) is 
defined as 



-1 



vn=0 



1 + /i 1 + /i a + 6 a 
1 — jj? 1 — 2ab I) 



(15) 



which agrees with more general expressions of Refs. 4, lo] 



III. SPECTRAL BIPHOTON WAVE FUNCTION 



Let us remind now the main features of the biphoton spectral wave function for SPDC 
with the type-I degenerate colhnear phase matching and with a pulsed pump. The main 
general expression has the form 

m 



where r is the pump-pulse duration, L is the length of a crystal, z/i and 1^2 are deviations of 
frequencies of the signal and idler photons cJi 2 from the central frequencies ujf'^ = 00^^ = 
ujq/2, |z^i,2| ^ ^0 is the central frequency of the pump spectrum, A and B are the 
temporal walk-off and dispersion constants 



Vg and Yg are the group velocities of the pump and ordinary waves, and ki and kp are the 
wave vectors of signal and pump photons. 

The first term on the right-hand side of Equation (|T6l) is the pump spectral field strength 
taken in the Gaussian form. The second term is the sinc-function characterizing the crystal 
and propagation in it of the pump and signle/idler photons. The argument of the sinc- 
function contains both linear and quadratic terms in 1/12- As we assume that 1 2^1,2! ^ 
in principle, the quadratic term is much smaller than the linear one. Nevertheless, the 
quadratic term can never be dropped because it is crucially important for determining the 
finite-width single particle spectrum [without the quadratic term both factors on the right- 
hand side of Equation f|T6l) depend only on z/i + i/2, which makes the integral / dh'2\'^{^i, ^^2)!^ 
independent of ui and single-particle spectrum infinitely wide]. As for the linear term in the 
argument of the sinc-function it can be efficiently eliminated in the case of long pump pulses 
but cannot be dropped if pulses are short. The control parameter separating the regions of 
short and long pulses is given by ^ 




(16) 




(17) 
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i.e., it's equal to the ratio of the double pump-pulse duration to the difference of times 
required for the the pump and idler/signal photons for traversing all the crystal. Pump 
pulses are short if <^ 1 and long if r/ 1 and, typically, ~ 1 at r ~ Ips. 

The degree of entanglement of the state ( fT6l) was characterized by two parameters, the 
Schmidt number K and the parameter R defined as the ratio of the single- to coincidence 
spectral widths of the corresponding photon distributions, R = ^u^'' / ^u^-f' [2]. Both pa- 

t[ps] 
2 3 4 5 




2 3 



n n 

FIG. 1: Entanglement parameters: R, calculated analytically [2], and K, calculated numerically 7], for a 
crystal LilOs of the length L — 0.5 cm and the type-I phase matching. 



rameters were calculated analytically as function of the pump-pulse duration r and were 
found to be rather close to each other. Also, for the same wave function ( |T6i) . the Schmidt 
number K{t) was calculated numerically [3], and the result was found to be in a very good 
agreement with the analytical function R{t) of Ref. ^ (see Fig. [1]). The degree of spectral 
entanglement was found to be rather high at all values of the pump-pulse duration and 
especially high in the cases of sufficiently short and long pulses (far from the minimum 
which was found to occur at ?7 ~ 1 or r ~ Ips). Schmidt modes were not calculated for 
the wave function of the form (JT6l) . either analj^ically or numerically. Below we present an 
approximate analytical derivation of the Schmidt modes separately in the cases of long and 
short pump pulses. 

For further analysis we will need explicit analytical expressions for the coincidence and 
single-particle spectral widths of the SPDC biphoton state f|T6l) in the limit of short pump 
pulses, 1] <^ 1. They were found in the work 2|] to be given by 



shotrt 



5.56 c 



is) 

shotrt 



2A\n{2)uJo 



(19) 
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Note finally that a very small difference between the Schmidt number K{t) (exact and 
numerical, and the parameter R{t) (analytical, [2]) raises often a question why are they 
so close to each other whereas the wave function (|T6l) is rather significantly non-double- 
Gaussian? In this work we give an answer: actually, the wave function ( |T6l) is "hiddenly 
double- Gaussian". And we make this hidden double- Gaussian character of the wave function 
(fT6|) explicit by finding the double-Gaussian models for the regions of short and long pump 
pulses and by generalizing these findings for all values of the pump pulse duration (Section 
6). 



IV. GAUSSIAN MODELING AND SCHMIDT MODES IN THE CASE OF 

LONG PUMP PULSES 



In the case of long pump pulses, rj ^ 1, the Gaussian function on the right-hand side 
of Equation ( JT6l) is much narrower the sinc-function. For this reason Ui ^ —1^2, and this 
eliminates efficiently the term linear in ui and z/2 in the argument of the sinc-function and 
reduces ^'(i^i, 1^2) to the form 

BL 



\E'(z/i, z/2) oc exp 



81n2 



smc 



2cuo 



(z/i - z/2)' 



(20) 



The wave function of a similar form was considered earlier in the analysis of angular en- 
tanglement of SPDC biphoton states jSj, and then some first Schmidt modes were found 
numerically. On the other hand, the form of the sinc-function in Equation ( l20l) is already 
very convenient for Gaussian modeling. Parameters of the corresponding Gaussian function 
can be chosen from the coincidence condition of the Full Widths at Half Maxima (FWHM) of 
the squared functions sinc(u^) and exp(— 7^^) which gives 7 = 0.249. Thus, the substitution 



smc 



BL 

2cujq 



Z/l - z/2 



exp 



-0.249 



BL 

2c Uq 



,1^1 - 1^2 



(21) 



reduces the wave function of Equation fl20l) to the double- Gaussian form ([T]) 

^2" 



^'(z/i, z/2) 



Tcab 



exp 



with the parameters a and b given by 



20ii2 



and b 



(z/i + Z/2)V^ 

2a2 



0.249 EL 



exp 



262 



27r 

0.249 ELAo 



(22) 



(23) 
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where Aq is the central wavelength of the pump. Note, that owing to the assumption 
?7 ^ 1 and also because always L ^ Aq, the ratio of the parameters h and a is large: 
b/a Tj^J Lj Ao ^ 1. By using Equations ([3]) we find also other parameters of the double- 
Gaussian wave function, /i and a: 



ong 



2 

CT 



0.249 In 2 SLA 







27r 



and 



ong 




r f0.2A9BL\o 
c \ 27rln2 



1/4 



The Schmidt number (ITSll for long pump pulses appears to be given by 



CT 



long 



1.5084- 



cr 



(24) 



(25) 



(26) 



2ab 2a 4^0.249 In 2 V^oLB V^oLB ' 

This result coincides perfectly with the earlier derived expression for the parameter R in the 



approximation of long pump pulses 



Vl-39 71 CT 



1.50739- 



cr 



(27) 



21n2 y/\;LB y/X^' 
Though the coefficients in Equations (l26l) and fl271) look different, numerically they are seen 
;o be amazingly close to each other. Note that the parameter -Riong was found in the paper 
2I without any Gaussian modeling and in the frame of a procedure absolutely different from 
that used above for the derivation of the Schmidt number -R'long- 

In terms of -ft'iong, Equations ( l24l) and ( l25l) for the parameters /x and a can be rewritten 

as 



1 



+ and «iong ~ = -7=L==. (28) 

^ long a a/ A long 2 a/ A long In Z 

Again in terms of -R'long, the Schmidt modes ffT^ and eigenvalues of the reduced density 
matrix (flU) take the form 



ong 



(i^iong In 2)1/4 



2WKiong In 2 



and 



A ^ 



lone 



1 



1 



long 



A. 



long 



long 



(29) 



(30) 



The derivations of this section were based on a direct Gaussian modeling of the wave 
function, without explicit consideration of the density matrix. As we will see in the following 
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section, in the case of short pump pulses such simple procedure does not work, and to apply 
Gaussian modeling, one has to calculate first the reduced density matrix. This indicates a 
rather important qualitative difference between the cases of long and short pump pulses. 

V. SCHMIDT MODES IN THE CASE OF SHORT PUMP PULSES 

In the case of short pump pulses we cannot use anymore the above-described procedure 
of a direct Gaussian modeling of the wave function ( JT6l) . Indeed, in this case the pump 
spectral function is much wider than the sinc-function and, hence, the linear term in the 
sinc-function's argument cannot be eliminated. As for the quadratic term, it cannot be 
dropped too because it is crucially important for appropriate determination of the single- 
particle spectrum. As for the sinc-function of Equation f|T6|) with its full argument, it's 
rather difficult to see directly any possibility of its replacement by any Gaussian function. 
For example, a simple substitution of sinc('u) by e~"^ gives a super-Gaussian function with 
quadratic and both third- and forth-power terms in 1/1^2 under the symbol of exponent. 
This is not a simplification at all and does not provide a way for analytical derivation of 
the Schmidt modes. So, we will apply here a somewhat subtler method. At first, we will 
calculate approximately the reduced density matrix corresponding to the wave function (fT6l) . 
The reduced density matrix found in such a way appears to be much more appropriate for 
double- Gaussian modeling than the wave function. For the Gaussian model of the reduced 
density matrix, its eigenfunctions (Schmidt modes), eigenvalues and the Schmidt number 
are easily found. After this, with the help of equations of Section II, we will find parameters 
of the effective double-Gaussian model for the original wave function (fT6|l . which appears to 
be very non-trivial and hardly guessable in advance. 
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A. Density matrix 



In a general form the density matrix corresponding to the wave function (fT6l) is given by 

[(z/i + Z/2)2 + {ul + Z/2) 



\2] ^2 



xsmc < — 



xsmc 



L 
2c 
L 



A{vi + z/2) - B 



A{u[ + U2) - B- 



81n2 



(31) 



Let us make in this equation a series of approximations similar to those made in Section 
VA of the work [2] and based, actually, on a single physical reason: under the assumption 
about short duration of pump pulses, r/ <^ 1, the Gaussian factor in Equation fl3T]) . as a 
function of either i^i, i//, or 1/2, is much smoother than both sinc-functions. Compared to the 
localization region of the Gaussian function ~ 1/r both the difference \vi — v{\ and sums 
|'^2 + 1^2 + ^/1 are very small: they are on the order of c/AL, as they are determined by 
the linear terms in the arguments of the sinc-functions. Owing to this, we can approximate 
z/2 by — z/i and by —v{ in small quadratic terms in the arguments of the first and second 
sinc-functions in Equation (!3T|) . thus linearizing their arguments with respect to V2'- 

[(z/i + v^f -f (z// + z/2)2] 



Pr{i^i, Vi) OC / dz/2exp 
AL 



81n2 



where 



xsmc 



l'2[l^l) 



2c 



[z/2 - z^2(z^i)]| sine N - '^2(«^l)]| 



_^45z/2 ,^4i?z//^ 
-z/i H — ^ and z/2(z/J = -z/^ + 



(32) 



(33) 



AijJq AijJq 

The slowly changing exponential function in Equation ( !32l) can be taken out of the integral. 



for example, at z/2 = | [z/2(z^i) + 1^2(^1)] to give 



exp 



r 



In 2 



,'\2 



16 



^^(z/i^ + z// 



/2^2 



exp 



r^5^(z/i + z// 
AA'^u'^ In 2 



/\4 



(34) 



where in the last approximate expression we have dropped all small terms proportional to 
(z/i — z//)^. In this approximation, instead of z/2 = | [i'2{^i) + '^2('^i)]) we could choose for 
evaluation of the exponential factor either z/2 = ^^2(^1) or z/2 = z72(z/]'). In all cases the final 
result would be the same as given by the last expression in Equation (1341) . 
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The remaining integral o 
Equation (37) of the paper 

di'2 sine 

AL 



the product of two sine-functions is calculated with the help of 



AL 



TT smc 



2c 



2c 



[u2 - z/2(z/i)] > sine 



AL 



'V2 - V2{vi)] 



TT smc 



AL 



(35) 



By combining Equations (IMll and ( l35l) together, we find the following representation for the 
reduced density matrix 



Pryi'i, I'l) oc exp 



4A2cugln2 



/\4 



smc 



AL 



[^1 



(36) 



Note that in the last approximate expression of Equation fl35|) all quadratic terms in the 
argument of the sinc-function are dropped. Though such approximation could not be used 
for the wave function f|T6|) . it appears to be well applicable for the density matrix fl35|) because 
the quadratic terms in the argument of the sine function are small and, most important, 
because the linear term is proportional to the difference of frequencies vi — v[ rather than 
their sum, whereas the exponential factor in Equation (1361) depends on vi + v{. This feature 
of the density matrix contrasts to that of the wave function (!T6|) in which both the pump 
spectrum and the linear term in the argument of the sinc-function depend on the sum of 
frequencies z/i + z/2. 



B. Gaussian modeling and Schmidt modes for short pump pulses 



Both factors in the reduced density matrix (l36l) can be modeled by Gaussian functions 
by means of substitutions 



and 



exp 



smc 



AA^ul In 2 



/\4 



exp 



-7i 



tB{ui + u[ 



i\2 



AL 
~2E 



exp 



-72 



2Au;oVhi2 
A^L\ui - u[f' 



(37) 



(38) 



where 71 and 72 are fitting factors to be defined below. 

With these substitutions the reduced density matrix (!36|) takes the double-Gaussian form 
® _ 



Pr[vi, z/J = - \ - exp 
a V vr 



15? 



exp 



262 



(39) 
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with 



In terms of the control parameter t] (fTSl) . roughly, 



«~^/^»^ (41) 

because for short pulses r] <^ 1 and always L ^ Xq. 

With the help of Equations ( |T3l) . dHD, and ( |T5i) we can find now other parameters of the 
double-Gaussian reduced density matrix (5 and /I), as well as the Schmidt number K. The 
latter is given by 



= ^ = (ln2)V^^ = ^-^(ln2)V^,/^. (42) 

C. Gaussian modeling of the original wave function 

In accordance with the discussion of Section 2, there is a one-to-one correspondence be- 
tween the double-Gaussian reduced density matrix and the double-Gaussian wave function. 
This means that, as we have received the double-Gaussian representation (!39l) for the re- 
duced density matrix of the biphoton state arising in the degenerate coUinear SPDC process 
with the type-I phase matching in the limit of short pump pulses, the original wave function 
( IT6|) also must be representable in the double-Gaussian form, even though it's not evident in 
advance. Parameters of the double- Gaussian model of the wave function can be found from 
Equations (HOj) and (jHj) , the latter of which can be simplified in the approximation a ^ 6 to 
give 



6^v^a = (ln2)i/%/^, (43) 



v/2 AL^ ' V ^1^^ 

with the relation between a and b opposite to that occurring between a and b (j4T 



a^b\j^<€.b. (44) 

With a and b given by Equations (H3|) . the double- Gaussian spectral wave function modeling 
that of Equation f|T6|) is given by Equation ([1]) (with xi^2 substituted by z/12). In the 
approximation a b (jHj) the coincidence and single-particle spectra corresponding to this 
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wave function are determined by the curves 



^1 



oc exp 



i2n 



V2=const 



and 



oc exp 



(45) 



(46) 



The coincidence and single-particle spectral widths are defined as FWHM of the curves (jl5 
and (iSl) 



2c /ln2 

AL' 



Az/f^ = bVh[2 = (In 2) 



3/4, 



'2Auo 



(47) 



72 ^ ' ■ y 7iT"^ 

By comparing these formulas with the earlier derived analytical expressions for spectral 
widths ( IT9|) . we find that all functional dependences in these two sets of formulas are identi- 
cal. As for numerical coefficients, they can be made identical too if we choose appropriately 
the fitting parameters 71 and 72 in the Gaussian model substitutes (1371) and ( l38l) for two 
factors in the reduced density matrix ( 136|) . ( l39l) . Found from these conditions parameters 
7i_2 are given by 

4 1n2 ln2 



71 



Vin2 ^ 0.832555, 



72 



0.0897. 



(48) 



(5.56)2 (2.78)2 

With these fitting parameters defined, we can write down now all final expressions for the 
model double-Gaussian wave function, Schmidt number, Schmidt modes and eigenvalues of 
the reduced density matrix for the case of short pump pulses. So, Equations (H3|) for the 
parameters a and b become equal to 

2.78 c 3.339 c , 



2AuJo 



.JWlAL AL ' y tB ' 
and the double- Gaussian model ([T]) of the wave function (1TB]) itself takes the form 



(49) 



' short 
72A2L2(z/i + Z/2^2 



AL 

TTC 



72 



K 



short 



X exp — exp 

2c^ 

0.045 A^L\ui + U2? 



7irE(z/i - z/2) 
AALUo^/\n2 



0.3 



short 



0.208 r5(z/i - z/2) 



X exp — '— exp \' . (50) 

Auj^^/Wl J ' 

Comparison of this expression with that of Equation (1T6|) shows that they ar not alike at all. 
But they do describe a very similar behavior of the wave function \E'(z/i, z/2). Similarity and 
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differences between the expressions ( !T6|) and (150|) are illustrated by two pictures of Figure 
|2j In these pictures localization regions of the wave functions f|T6l) and fl50l) are shown by 
thick solid lines. "Centers of mass" of these lines correspond to maxima (equal to unit) of 
the narrowest factors in the formulas for these wave functions. For the exact wave function 
( TT6l) the narrowest factor is the sinc-function, which is maximal when its argument equals 
zero. This condition yields 

^Jexact (^2) = ^2 + ^ (a - ^ + 9>BAu,/uJ,) . (51) 

^^shotrt ~ 5.56c/y4L. For the model double-Gaussian wave function the narrowest part is 




FIG. 2: Thick solid lines determine localization regions of the wave functions (a) (fT6|) and (6) (f50|) . dashed 
lines determine limitations of these regions along the thick solid lines; calculations for LilOg, pump wave- 
length Ao = 400 nm and pump-pulse duration t = 50 fs. 

the first exponential factor on the right-hand side of Equation (!50l) and, hence, 

= -^2. (52) 

Widths of the thick lines in Figure [2] are determined by the coincidence spectral width 
for short pump pulses ( fT9i) . /S.vf' = 5.56 c/AL, and, thus, the wave-function localization 
regions are regions under the thick solid curves in the (z/i, z/2)-map. Note, that in the 
pictures of Figure [2] widths of the thick solid curves are slightly increased compared to At/{ 
to distinguish clearer these curves from other lines in figures. 

Dashed lines in the pictures of Figure [2] determine limitations of the wave-function lo- 
calization regions in the longitudinal direction imposed by the remaining slowly varying 
factors in Equations ( JT6l) and ( |50l) [the Gaussian pump spectral function in ( JT6l) and the 
second Gaussian factor on the right-hand side of Equation (!50|) ]. Mathematically these lim- 
itations can be found from the condition that the corresponding squared Gaussian function 
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in Equations (fT6|) and (150|) are larger or equal to, e.g., 1/2, which gives 



cm. ^ dashed/,, \ , .dashed — z, , N ^ , ,c.m. \ ^ , .dashed + , N ^CQ^ 

1 exact (^^2) < '^l exact (^^2), '^l Gauss < '^l Gauss ( ^^2 ) < '^l Gauss , (53) 



where i^iexact (^2) and i^iGauss^('^2) are the functions determining location of the dashed lines, 
correspondingly, in Figures [2]^ a) and[2]^6): 



-S^(-2) = -.2 + «s^(-2) = .2 ± (ln2)^/y^;^. (54) 

These formulas, as well as the pictures of Figure [21 show clearly that mechanisms, limiting 
the localization regions of the exact f|T6|) and model fl50l) wave functions are absolutely 
different. But results are seen to be almost identical. Moreover, the only difference in the 
resulting localization regions in the pictures (a) and (6) of Figure [2] is a very small curvature 
of the thick solid line in Figure [2](a) missing in Figure [2t^6). This curvature arises owing to 
the quadratic term in the argument of the sinc-function in wave function f|T6|) . As discussed 
above, this term is crucially important for determining the single-particle spectrum and 
limitation of the localization region of the exact wave function (fT6!) . In the model function 
( 150|) quadratic terms are missing. Moreover, even if taken into account in some way, such 
terms are not expected to be important at all. All integral characteristics of the SPDC 
spectra are identical if calculated with the exact and model double-Gaussian wave functions. 
In particular, this is true for for the Schmidt number K and single-particle and coincidence 
spectral widths. As for the Schmidt modes and eigenvalues of the reduced density matrix, it 
would be very interesting to calculate them numerically to compare with the above derived 
formulas (l58l) and (1601) . We hope to be able making such comparison later. 

With the fitting parameters 71^ 2 determined by Equation ( l48i) and the parameters a and 
h given by Equations (H9|) . Equation fH2l) for the Schmidt number -ft'short takes the form 



^ ± ^ V^tal ^ i ^ 0.5308 . (55) 

2a 2.78 y^yO^ VB VXi^ ^ ' 

f' 

This expression coincides exactly with that of the work [2] for the parameter R in the regime 
of short pulses. 

The reduced density matrix fl39|) . after all substitutions, takes its final form 



Pr,short(,^lj ^1; 



X exp 



Br , 







exp 



2Bt 

^2^2 In 2, ,,2' 



(5.56 c) 



(56) 



16 



For the parameter /i, determining eigenvalues of the reduced matrix, we get formally the 
same expression as in the case of long pump pulses, ~ — 1 + 2|, though with different 
constants a and b, which gives finally 

^^i^short" ^^V^AV^ L ' ^^^^ 
and eigenvalues of the reduced density matrix are given by 

An-T^fl-T^V (58) 



-^short \ -^short , 

with -K'short determined by Equation fl55|) . 

The scaling factor a of the Schmidt modes (fT3|) in the case of short pump pulses has the 
form 

1 AL Vh[2 AL , , 

a.hort = = ^ = 0.299 , (59) 

a V -f^ short C 2.78 V-f^short CV-f^ short 

and the Schmidt modes are given by 



, , AL (In 2)1/4 vi,2AL Vln2 

^"^'^^'^^ - ^/ V7!:ff(W^''" i C 2.78Vi^short 



0.547,/^^ V^Jo.299^^). (60) 

CV-f^ short \ CV-f^ short. 



VI. GENERALIZATIONS 

By comparing Equations fl2^ and fl5Ul) with flFTl) and fISS]) . we find that they look identical 
with the only substitution: -ft'iong ^ -f^short- So, we can suggest the following simplest 
interpolation rule for a transition between the regions of short and long pulses consisting in 
the assumption that even in the intermediate region {rj ~ 1) the dependence of and A„ of 
the Schmidt number K remains the same as in the asymptotic regions rj and r/ ^ 1: 

/i(r) = -1+ ^ 



K{t) 



Of course, these relations are proved rigorously only in the limits of short (r^ <^ 1) and 
long [r] ^ 1) pump pulses. But as the Schmidt number K{t) is large also in the case of 
intermediate pulse durations, the assumption that Equations ( !6TI) remain valid also in the 
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case ?7 ~ 1 seems rather natural. Note however that the function K(t) in the region of 

n 

intermediate pulse durations either has to be taken from numerical calculations [7] or in its 
turn has to be defined by means of a more or less reasonable interpolation. The simplest 
interpolation K{t) = \/ K^i^orti'^) + ^^lonfri'^) suggested and used in our earlier work j^. 



long \ 

Below we discuss somewhat more elaborate ways of interpolating K{t) into the region ~ 1. 

Following the same logic, let us make now a rather important assumption that the double- 
Gaussian modeling of the wave function f|T6|) has sense not only in the asymptotic cases of 
short and long pump pulses but also in the case of pulses of intermediate duration, when 
?7 ~ 1. In other words, we assume that for any values of the control parameter t], small, 
intermediate, and large, the biphoton wave function (ITB]) can be satisfactorily modeled by 



^(1^1,1^2; r) = J — exp 
7ia[T)b[T) 



[V1 + V2? 



2a(r) 



exp 



2b{T) 



(62) 



where the widths a(r) and 6(r) are determined by Equations ( 123|) and ( l49i) in the cases of 
long and short pump pulses and have to be defined yet in the case of intermediately long 
pulses. By comparing expressions in Equations (l23l) and fH9l) we find easily that they are 
related to each other by the following simple formulas 

flshort 1-39 



= 1.00267 77 ~ and 



along 2 In 2 

1.00207 ^ ^. (63) 



^long _ / CT 



&short V 2 X 0.249 LA 

By using the idea of interpolation we assume that in a general case we can define the 
functions a(r) and a(r) as being given by 

air) = a^hovtfaiv) and 6(r) = konghiv), (64) 
where the smoothing functions obey the conditions 



/a(^)lr,«l = 1. /a(^)l„»l = -; /fe(^) = /b(^) = 1- (65) 



Of course, there are infinitely many ways of choosing the smoothing functions fa{v) and 
fbiv)- To restrict somehow this manifold of possibihties, let us take these functions in the 
form 

/„(r7) = (l + r7^)"^/^ and /.(r^) = ii±^!L^, (66) 
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where s > is a free parameter. With these functions faiv) hij]) Equations (167|) are 
reduced to 



(1 + „.)-"• ;^^ = -i:f=^(l + „r'". 



b{T) = 6long 



[1+7]' 



bir) 



(67) 



^ ' tuo V 27r 0.249 EL ^ 

The dependences a(r) and 6(r) are characterized by two curves in Figure [3l The dashed parts 
of these curves characterize smooth transitions from the short-pulse to long-pulse regions. 
One of the main features of the widths a(r) and 6(r) seen clearly in Figure [3] is that at all 
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FIG. 3: Parameters a(r) and 6(t) (IST)) of the general double- Gaussian spectral wave function calcu- 
lations for LilOs crystal, L — 0.5 cm; dashed lines are interpolation with the help of the functions (pSj) with 
s = 2.21. 

values of the pump-pulse duration r, or of the control parameter rj, always 6(r) ^ a(r). 
This feature of the model wave function ([H2D differs significantly from that of the exact 
wave function (ITB]) . In the latter the pump spectral function and sinc-function change their 
roles at ?7 ~ 1: in the region <^ 1 the sinc-function is much narrower and in the region 
1] ^ 1 much wider than the pump spectral function. In contrast to this, in the case of the 
model double-Gaussian wave function, the first factor on the right-hand side of Equation 
f l62|) is much narrower than the second one at all values of the pump-pulse duration r. This 
simplifies, for example, the definition of the coincidence and single-particle spectral widths. 
In terms of a(r) and 6(r) they are given by 



Ai/J'^ (r) = 2a(r) 7^2, Au['^ (r) = 6(r) 7^2 



(68) 



for all r. 
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With the definitions of Equations (167|) we get the following generalized expression for the 
Schmidt number 

K(r) = ^= + 

^ ' 2a(r) 2.78V2 x 0.249 V BX^ ^ ^ 

The dependence K{t) (or K{r])) is similar to that characterized by the curves in Figure [H 
The minimum of the function K(t) is achieved at r]o = 2^^/* and 

^^^^'-''^ " 2.78V 2 X 0.249 ^ V 5^ 2^" ^^"^ 
The case s = 3 corresponds to the earlier used interpolation j^. If we want to increase 
^i'^)\mm -'-2^ make the curve K{t) closer to the numerically calculated one 
we have to take s = 2.21 in our formulas (lM]) - (|7n]) . This gives K{t)\^^^ = 83 (instead 
of K{t)\^^^ = 73 of the work and the curve K{ri) appears in this case practically 
indistinguishable from the numerically calculated one. But of course, the main novelty 
of the described here method of interpolation is not only in a possibility of getting this 
12% improvement. The main message is in the assumption about a possibility of modeling 
the wave function (fT6|) at arbitrary values of the pump-pulse duration r by the double- 
Gaussian wave function fl52]) with the parameters a(r) and 6(r) given by Equations fl^TI) . 
This assumption is proved to be correct in the cases of long and short pulses. At the 
"integral level", correctness of the suggested model (l62l) is also proved for all values of r 
because the Schmidt number K{t) is found to be practically identical in the exact numerical 
calculations ^] and in the given above analytical derivation for the model double-Gaussian 
wave function fl62|) with s = 2.21. Localization regions of the exact f|T6|) and model fl62|) 
wave functions in the case of intermediate pump-pulse durations are also very close to each 
other and, qualitatively, similar to that shown in Figure [2](6). A more detailed comparison of 
the exact f|T6l) and model fl62l) wave functions in the case of intermediately long pump pulses 
can involve comparison of eigenvalues of the reduced density matrix and of a structure of 
the Schmidt modes. But for the exact wave function f|T6l) at r/ ~ 1 such numerical data are 
not obtained yet. 

Returning to the generalized formula (IMl) for the eigenvalues A„(r) of the reduced density 
matrix, as said above, we assume that this expression is valid for not only for small and 
large t] but also also in all the region between these two asymptotic limits. In accordance 
with the Schmidt theorem, the weights with which the Schmidt modes are represented in the 
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expansion (fTTj) are given by ^/KAj) 
K{t) = 87. In this case the values of 
H At n = K = 87, we get v^VF) 



. As an example, let us consider the case rj = 1, when 
a/ A„(r) are located along the hne shown in Figure 



r?=l 



r,=l 



0.055, which as high as 37% of the maximal value 
of |-\/A^|, = 0.15. This results show that falling of |-\/A^| with a growing n is rather 

slow, and in reality there is a rather large interval of n > , giving not too small contribution 
into the expansion of the wave function in the sum of the Schmidt-mode products (JTTj) . 
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FIG. 4: Square roots of the eigenvalues of the reduced density matrix (|6ip at 77 = 1 and K{t) given by 
Equation 



The Schmidt modes of the model double-Gaussian wave function fl62l) have the same 
form as in the asymptotic cases of long (129|) and short (!60|) pump pulses, but now with the 
generalized parameters a(r) and 6(r) ( 1671) : 



1 



a(r) 



uJoa[T) 



^1,2 



(71) 



where u;oQ;(r) is the dimensionless generalized scaling factor characterizing the dependence 
of ipn on the dimensionless variables z/i^ 2/1^0: 

3/4 



a;oa(r) = ujq^ 



/ r \ 3/4 



a{T)b{r) 

As a function of the control parameter 77, the dimensionless scaling factor of the Schmidt 
modes is plotted in Figure [51 The dimensionless scaling factor is large practically at all 
values of rj. With a growing rj the scaling factor monotonously grows, which means that 
localization regions of all Schmidt modes shrink. At rj = 1, Eq. ( IVII) gives ujQa{T) ^ 585.5. 
The spectral width of an n^^ Schmidt mode (17T1) can be estimated as 



n 



a(r) 



n a(r)6(r) 



(72) 
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FIG. 5: Dimensionless scaling factor woc«(t) (jVip vs. the control parameter 77, for LilOa crystal, L = 0.5 cm, 
Ao = 400 nm, s ^ 2.21. 



This width grows with the increasing mode number n as ^/n. For n = K{t) Equation (ff2ll 
gives 6uk ~ &(r)/2 = Az/(^V2Vin2 = 0.6 Az/(^); 5z/„(r) reaches Az/(^) at n ^ 2.77 K{t). It's 
true, however, that the contribution of such high-number Schmidt modes to the expansion 



(fTTj) is very small: for the same parameters as used above for all estimates, we get \f\2riK ~ 
6% of V%" (to be compared with ~ 37% of \f\\). This means that practically for all 
important Schmidt modes their spectral width is smaller than the single-particle width of 
the biphoton spectrum as a whole. 

As an example, some spectral Schmidt modes are shown in Figure E] for the same param- 
eters as in Figure O and r/ = 1. 
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FIG. 6: Schmidt modes V'lj V'lo and t\)' 
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VII. CONCLUSION 



Summarizing, we found double-Gaussian models of the SPDC type-I wave function fll6p 
in both regions of short and long pump pulses. In the case of short pump pulses even 
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a possibility of Gaussian modeling was not evident in advance. Its derivation required a 
transition to the reduced density matrix where Gaussian modeling appears to be much more 
natural and simple than in the case of the wave function. Then, after finding a double- 
Gaussian model for the reduced density matrix and using general relations between its 
parameters and those of the wave function, we were able to find the double-Gaussian model 
( !50|) of the original wave function (fT6l) in the case of short pump pulses. Though the forms 
of these too functions are significantly different, the double-Gaussian model reproduces a 
structure of the exact wave function rather well and almost in all details, which is seen 
clearly in two pictures of Figure O Having the double-Gaussian models of the exact wave 
function in the regions of short and long pulses we found also in these two cases the Schmidt 
modes ( 129|) and ( l60l) . Finally, by analyzing features of all parameters in the asymptotic 
cases of long and short pump pulses and using the most appropriate interpolation to the 
region of intermediately long pump pulses, we found a possibility of suggesting the model of 
a double-Gaussian wave function [Equations ( !62l) and ( !67l) ]. which we assume to be valid for 
arbitrary values of the pump-pulse duration r. For this model wave function we found also 
the Schmidt number K{t) ( l69l) . Schmidt modes [Equations ( 17T1) and (IVII) ]. and eigenvalues of 
the reduced density matrix fl6T|) . The Schmidt number found in such a way agrees perfectly 
well with the results of its numerical calculation We assume that the functions ipn{i^i,2) 
of Equation fl7T|) represent equally well the spectral Schmidt modes in all range of the pump- 
pulse durations, short, long and intermediate. Comparison with exact numerical calculations 
of the Schmidt modes is expected to be interesting and fruitful. Also we hope that the derived 
results can be tested experimentally. For example, an interesting scheme of an experiment 
can involve splitting of the biphoton beam for two channels by a nonselective beam-splitter, 
installing in one of two channels a spectral mask gating through it only a single spectral 
mode and measuring the coincidence spectrum in the second channel. According to the 
Schmidt theorem the coincidence spectrum measured in such a way is identical to that of 
the mode gated through in the first channel, and this coincidence of two spectra deserves 
its experimental verification. We hope also that in future the data about spectral Schmidt 
modes can be used for practical applications in the problems of quantum information and 
quantum cryptography. 
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